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We have examined the production of J/tp mesons in high energy proton- proton and proton- nucleus 
collisions at beam energies in the range from 158 GeV to 920 GeV, available from different fixed 
target experiments. In the employed model J/ip production in hadronic collisions is assumed to 
be a factorisable two step process: (i) production of a cc pair which can be reliably described by 
perturbative QCD, and (ii) formation of J /if) resonance from the cc pair, which can be conveniently 
parameterized incorporating different existing physical mechanisms of color neutralization. We 
show that, for lower collision energies, J/ip production through quark-anti-quark annihilation gives 
larger contribution at higher xf, while gluon-gluon fusion dominates the production at smaller 
Xf- For proton-nucleus collisions the model takes into account both the initial state modification 
of parton distributions in nuclei and the final state interaction of the produced cc pairs with the 
target nucleons. The model is found to give reasonable description of data on J/ip production 
in protonic and proton-nucleus collisions, for different existing fixed target experiments. In case of 
proton-nucleus collisions, our calculations show a non-negligible dependence of the final state nuclear 
dissociation of J/ip mesons on the energy of the incident proton beam. The model has been applied to 
predict the J /if) production and suppression expected in proton-nucleus collisions at energies relevant 
to FAIR, the upcoming accelerator facility at Darmstadt, Germany. The amount of suppressions, 
for different mechanisms of 3/ip hadronization has been found to be distinguishably different which 
might help an experimental settlement of the much controversial issue of color neutralization. 



I. INTRODUCTION 



Collision of heavy nuclei at relativistic energies endows 
us with the unique opportunity to create and investigate 
hot and dense nuclear matter in the laboratory. Statisti- 
cal calculations of quantum chromodynamics (QCD), the 
fundamental theory governing the dynamics of strongly 
interacting particles predict that at sufficiently high tem- 
perature and/or density nuclear matter undergoes a tran- 
sition from the confined hadronic phase to a color de- 
confined partonic phase. Hence the primary objective 
of the relativistic heavy-ion physics program is the ex- 
ploration of the phase diagram of the strongly interact- 
ing matter, in various regions of temperature and baryon 
density and to confirm the existence of this novel phase 
of quark matter, the quark-gluon plasma (QGP). Over 
the last three decades relentless efforts have thus been 
invested to devise unambiguous and experimentally vi- 
able signals that would clearly identify the existence of 
this phase transition and the signatures of this tran- 
sient phase. The suppression of charmonia in high en- 
ergy heavy-ion collisions, has long been predicted as a 
'smoking gun' signature for the de-confinement transi- 
tion PHI- More than 25 years ago, Matsui and Satz ar- 
gued in their seminal paper, that the production yield of 
J/ip mesons, a bound state made up of charm anti-charm 
pair, should be considerably suppressed by Debye screen- 
ing of the colored partons if relativistic collisions would 
result in the formation of de-confined QCD medium [l|. 
However subsequent experimental investigations have re- 
vealed a considerable reduction of the charmonium pro- 
duction already present in proton-nucleus (p-A) colli- 
sions compared to the hadronic collisions (scaled with 
binary nucleon-nucleon (NN) collisions) [3t13 • A com- 



plete understanding of this 'normal nuclear suppression' 
is mandatory to establish a robust unambiguous baseline 
reference , with respect to which the 'anomalous suppres- 
sion', pertaining to formation of a de-confined partonic 
medium can be clearly identified. 

The upcoming compressed baryonic matter (CBM) ex- 
periment at the Facility for Anti-proton and Ion Research 
(FAIR) aims at the exploration of the QCD phase dia- 
gram in the region of high baryon densities and mod- 
erate temperatures, by colliding light and heavy-ions in 
the beam energy range E Lab = 10-40 AGeV HJ. Model 
calculations performed so far [l4| , based on transport as 
well as hydrodynamical equations, indicate that for cen- 
tral (b=0) collisions highest net baryon densities ps = 6 
- 12 pq (where po is the density of normal nuclear matter) 
are expected to be produced at the center of the collision 
zone. The experiment will thus provide us the unique 
possibility to create and investigate extremely dense nu- 
clear matter in the laboratory through measurements of 
bulk as well as rare probes. Such measurements will en- 
able us to have deeper insight on some of the fundamen- 
tal but yet enigmatic issues of strong interaction ther- 
modynamics like, hadronic properties in dense nuclear 
matter, de-confincmcnt phase transition from hadronic 
to quark-gluon matter driven by high baryon densities 
and the nuclear equation of state (EOS) at high baryon 
densities. The charmonium produced at an early stage 
of the collision is believed to be a penetrating probe for 
characterization of the confining status of the dense bary- 
onic medium. The experimental program thus includes a 
detailed proposal to study the production of the hidden 
charm hadrons, through their decay into lepton pairs [l5j . 
by high intensity beams. The di-leptons are not affected 
by the strong interactions that reign during hadroniza- 
tion and carry undistorted images of the earlier phases. 
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At FAIR, production of charmed particles will be stud- 
ied at energies close to their kinematic threshold (the 
threshold beam energy for 3/tp production in p-p colli- 
sions being E t h ~ 11.3 GeV) and hitherto unexplored by 
the existing experiments, due to extremely low produc- 
tion cross sections. 

In this paper we present a detailed study of the char- 
monium production in proton-proton (p-p) as well as 
proton-nucleus (p-A) collisions and its observed suppres- 
sion in case of nuclear targets. To account for the nuclear 
reduction of the 3/tp yield, some of the so called 'cold nu- 
clear matter' (CNM) effects affecting the production in 
p-A collisions, have been incorporated in the model. Re- 
sults from our calculations have been compared and con- 
trasted with the data for 3/tp production at fixed target 
energies. The outcomes of our investigation have been 
utilized to predict the 3/tp production and its suppres- 
sion expected to occur in the FAIR energy domain. 

The paper is organized in the following way. The model 
for calculation of 3/tp production in hadronic as well as 
nuclear collisions is briefly depicted in section II. Section 
III describes the analysis of the existing 3/tp production 
cross section data in p-p and p-A collisions, available 
from different fixed target experiments. In section IV 
we present our results for the estimation of 3 ftp yield ex- 
pected in energetic collisions at FAIR. Finally in section 
V we summarize our results. 



II. THEORETICAL FRAMEWORK FOR 
CALCULATION OF 3 /ip PRODUCTION 

In hadronic collisions, one parton (quark or gluon) 
from each of the colliding hadrons interacts strongly to 
produce a cc pair. Subsequently the cc pair hadronizes to 
form a physical 3/ip meson, provided the invariant mass 
of the pair is below the threshold for open charm produc- 
tion. Our present work has adopted the model originally 
proposed by Qiu, Vary and Zhang [1 6j j and subsequently 
used by others [17| for interpretation oi 3/tp suppression 
in nucleon-nucleus and nucleus-nucleus collisions at SPS 
and RHIC. In this model, the 3/tp production in high 
energy hadronic collisions is believed to a factorisable 
two stage process. The first stage is the production of a 
cc pair with relative momentum square q 2 . Because of 
large charm quark mass, the process involves a short time 
scale and can be well accounted by perturbative QCD. 
The second stage involves the formation of the color neu- 
tral physical bound state 3/tp, from the initially compact 
cc pair. It takes a relatively longer time and thus non- 
perturbative in nature. Over the years, the main debate 
on the mechanism of 3/ip production is focused on this 
second stage, which is not yet fully understood on a fun- 
damental theoretical level. Three models are commonly 
found in literature for calculating the cross sections of 
3/tp production: the Color Evaporation Model (CEM) 
PJ, the Color Octet Model (COM) and the Color 
Singlet Model (CSM) [2(J. The production rate of the cc 



pairs with an invariant mass Q 2 , can be factorized into 
i) a convolution of two parton distributions from the two 
incoming hadrons and, ii) da a +b^>cs+x / 'dQ 2 , which rep- 
resents the perturbatively calculable short-distance hard 
parts for the parton a and b to produce the cc pairs with 
mass Q 2 [2l[. At the leading order, the partonic contri- 
butions come from two subprocesses: quark annihilation 
(qq — > cc) and gluon fusion (gg — > cc). With the K- factor 
accounting for effective higher order contributions, the 
single differential 3/ip production cross section in colli- 
sions of hadrons h\ and ft. 2 , at the center of mass energy 
\fs can be expressed as, 



da 



j/il> 

hi hi 



dxi 



= K 



dQ 2 



da' 



h\h<2 



dQ 2 dx J 



;) x F^j/^q 2 ), (1) 



where Q 2 = q 2 + 4m^ with mc being the mass of the 
charm quark and xf is the Feynman scaling variable. It 
can be related to x a and Xb by xf = x a — Xb, with x a and 
Xb being the momentum fractions carried by the incom- 
ing partons originating from the beam and the target re- 
spectively. In Eq. (TTJ F cS ^j^(q 2 ) is the transition prob- 
ability that a cc pair with relative momentum square q 2 
evolve into a physical J/ip meson, in hadronic collisions. 
Three different parametric forms have been formulated 
for the transition probability corresponding to the three 
above mentioned physical mechanisms of color neutral- 
ization: 
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where mo is the mass scale for the open charm produc- 
tion threshold (DD threshold). In F cS ^j/^(q 2 ), Nj/^ 
and qf are the model parameters that can be fixed by 
comparing the model results with the existing total pro- 
duction cross section data from hadron-hadron collisions. 

As illustrated in [lj|, the different functional forms of 
the transition probability represent different pictures of 
color neutralization. The F^(q 2 ) implies that a con- 
stant fraction of all the produced cc pairs with invariant 
mass below the open charm production threshold evolves 
into the physical 3/tp mesons, and thus bears the central 
theme of the Color Evaporation Model [U[ . F^ G ) (q 2 ) car- 
ries the essential ingredients of the Color-Singlet Model 
[20| . It is assumed that the cc pair is produced in a 
color-singlet state. The transition amplitude (cc\3/ip) 
does not involve any radiation and interaction with the 
medium, and thus proportional to square of the 3/ip wave 
function parameterized as Gaussian. Finally, F^(q 2 ) in 
Eq. Q, mimics the the essential features of Color-Octet 
Model [19]. In this model the cc pairs are assumed to 
be produced in color-octet stage. Subsequent formation 
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FIG. 1: The MSTW 2008 leading order (LO) parton distribution 
functions in a free proton for gluon (continuous line), up quark 
(long dashed), down quark (short dashed), anti-up quark (dotted 
line) and anti-down qurak (dashed-dot line). The parton densities 
are evaluated at a scale Q = Wlj/y,. 



of color-singlet physical resonances occurs through the 
gradual expansion of the initially compact cc pairs. The 
expansion is associated with the radiation of soft gluons 
necessary for color neutralization. The g 2 -dependence 
of the transition probability is assumed to be associated 
with that radiation, and a power-law (P) distribution, is 
believed to represent the transition probability. 

Instead of Xf, one can also express the single differen- 
tial cross section, in terms of the center of mass rapidity 
Vcms-, of the cc pair as, 



da 



J/4> 

hih 2 



(X a + X b )- 



da 



j/i> 

h±h 2 



dxF 



(5) 



The inclusive cross sections as reported by different 
experiments can then be estimated by integrating over 
corresponding kinematic range. 

The double differential cross section for producing a cc 
pair can be decomposed into the individual contributions 
from different channels and at leading order is given by 



da 



h\h 2 



do 



h\h 2 



da" 
aa h 1 h 2 



dQ 2 dxF dQ 2 dxF dQ 2 dxF dQ 2 dxF 



(6) 



where a qq is the contribution from annihilation of light 
quarks, q = u,d, s, a CCint is the contribution from intrin- 
sic charm quarks, and a" is the contribution from fusion 
of two gluons. Due to very small charmed sea distribu- 
tion appearing quadratically in the cross section calcu- 
lation, the contribution of intrinsic component a CCint , is 
negligible compared to the other terms |22j. Hence the 
significant contributions from the other two leading terms 
are 



for light quark annihilation, and 



da" 

aa h x h 2 

dQ 2 dx} 



T"{Q 2 ) XgX b 

Q 2 ; 



Xb 



ghi(x a ,Q 2 )g h2 (x b ,Q 2 ), 



(8) 

for gluon fusion. Here q(x,Q 2 ), q(x,Q 2 ) and g(x,Q 2 ) 
signify the quark, anti-quark and gluon distributions, re- 
spectively, in a hadron, evaluated at a scale Q 2 . In the 
present work, we have used MSTW 2008 [I| LO par- 
ton distribution functions for free protons. Fig. Q] shows 
the momentum distributions for up and down quarks, 
anti-quarks and gluons following MSTW 2008 LO central 



set, evaluated at the scale Q 



in 



Because of the 



two-parton final-state at the leading order, the incom- 
ing parton momentum fractions are fixed by the kine- 
matics, and at a given i/s, can be expressed in terms 
of Xf and Q 2 , 



by x a = (\/x 2 F + AQ 2 /s + xf)/2 and 
xp)/2, respectively. The partonic 
cross sections appearing in Eqs. (J7]) and flHJ, in the lead- 
ing log approximation, are given by (22l . [23( 



.n, = {^fx 2 F + AQ 2 /s 
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and 



(Q 2 ) = gf[(l+7 + ^7 2 )K 



/ 7 31 r— 



7. 



(10) 



where a s is the QCD running coupling constant, and 
7 = Am 2 c /Q 2 . 

Let us now turn to a description of 2 production in 
proton-nucleus collisions. In this case, the presence of 
normal nuclear matter can affect the charmonium pro- 
duction and thus these collisions provide a tool to probe 
the effect of confined matter. Nuclear effects can come 
into play during the entire evolution period of pro- 
duction. Several different phenomena collectively known 
as cold nuclear matter (CNM) effects have been studied 
in literature at considerable detail. Two most important 
CNM effects that have been identified are: i) the ini- 
tial state nuclear modification of the parton distribution 
function affecting the perturbative cc pair production, ii) 
the final state dissociation of the nascent cc pair in the 
pre-resonance or resonance stage, due to its interactions 
with nucleons during its passage through the target nu- 
cleus. 

The PDF in a nucleus with atomic number Z and mass 
number A is written as the sum of the proton (ff ) and 

the neutron (f^ A ) parton densities in a nucleus: 



ft A =z ff /A +(A-z)f: /A , 



(11) 



dQ' 2 dx p- 



Z qq (Q 2 ) x a x b 
Cp x a +x b 



where is obtained from ff /Ji by isospin conjuga- 



p/A 



tion: it 



n/A _ ^p/A ^n/A _ U p/A g n/A 
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Deep 



x [qhi^ (x a ,Q 2 )qh2H x b,Q 2 ) + qiii? (x a ,Q 2 )qh 2 f ( x b,Q 2 Jft) inelastic scattering (DIS) and Drell-Yan measurements 
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FIG. 2: Nuclear modification of valence up quark (long dashed 
line), sea up quark (short dashed line) and gluon (continuous line) 
distribution inside a Pb nucleus following EPS09 interface evalu- 
ated at the scale Q = mji^. EPS09 gives Ri(Pb,x,Q 2 ), the scale 
dependent ratio of the distribution of the parton flavor i for a pro- 
ton in a nucleus A, to the corresponding parton distribution in the 
free proton. 

performed with nuclear targets have shown that the dis- 
tributions of partons in nuclei are significantly modified 
relative to those in free protons. These nuclear modifica- 
tions depend on the fraction of the total hadron momen- 
tum carried by the parton, x, on the momentum scale, 
Q 2 , and on the mass number of the nucleus, A. While 
the mechanisms governing these modifications are not yet 
well understood, several groups have produced param- 
eterizations, Ri(A, x, Q 2 ), that convert the free-proton 
distributions for each parton i, ff{x,Q 2 ), into nuclear 
ones, f^(x 7 Q 2 ), assuming factorization: 

tf{x,Q 2 ) = R l (A,x,Q 2 )* f?{*,Q 2 ) ■ (12) 

For our calculations, we have employed recently pro- 
posed EPS09 [111 interface, which is available for all mass 
numbers. Fig. [2] shows the ratio Ri(A,x,Q 2 ) following 
EPS09 parameterization, in a Pb nucleus for valence up 
quarks, sea up quarks and gluons, calculated at the scale, 
Q = mji^. It is important to note that, the quark and 
anti-quark distribution functions are directly probed by 
the nuclear DIS and Drell-Yan data, their nuclear effects 
are relatively well constrained and different parameteri- 
zations give almost similar results. The connection be- 
tween the measurements and the nuclear gluon densities 
is much more indirect. Since gluon fusion plays the dom- 
inant role for cc production, a complete understanding 
of charmonium production in proton-nucleus collisions 
is presently lacking due to uncertainties in nuclear gluon 
distributions, as illustrated in [26| . It may be mentioned 
here that in the original model by Qui, Vary and Zhang 
the nuclear effect to the parton distribution functions are 
ignored and free proton pdfs have been used to estimate 
the cc pair production in nucleon-nucleus and nucleus- 
nucleus collisions. 

For simulating the final state interactions we closely 
follow the method advocated by Qui, Vary and Zhang 



[la . |22( , who have treated the conventional nuclear dis- 
sociation in somewhat unconventional way. In nucleon- 
nucleus collisions, the produced cc pairs are likely to in- 
teract with the nuclear medium before they exit. Ob- 
served anomalous nuclear enhancement of the momen- 
tum imbalance in di-jet production [27j indicates that a 
colored parton (quark or gluon) experiences multiple soft 
scatterings when it exits the target nucleus (28[. Hence, 
the colored cc pair produced through partonic hard scat- 
tering should also experience multiple soft scatterings 
when it passes through nuclear matter. These multiple 
soft scatterings will increase the relative transverse mo- 
mentum between the c and c, and consequently, increase 
the invariant mass of the cc pair. Some pairs can thus 
gain enough momentum to be pushed over the threshold 
and become two open charm mesons, and the J/ip pro- 
duction cross section will in turn be reduced in compar- 
ison with nucleon targets. Larger be the size of the nu- 
cleus, more be the number of soft scatterings undergone 
by the pair and the reduction in the cross section for J/?/> 
production will be even greater. If the formation length 
for the J/ip meson, which depends on the momenta of the 
cc pairs produced in the hard collision, is longer than the 
size of the nuclear medium, it is reasonable to assume 
that the transition probability F cS ^j/^,(q 2 ), defined in 
Eq. ([T]), can be factorized from the multiple scattering. 
Then, as far as the total cross section is concerned, the 
net effect of the multiple scattering of the cc pairs can be 
represented by a shift of q 2 in the transition probability, 

q 2 — > q 2 = q 2 + S 2 L(A) . (13) 

In Eq. (fT3|). L(A) is the effective length of nuclear 
medium traversed by the cc pair, from its point of pro- 
duction till it exits and depends on the details of the 
nuclear density distributions (29|. In Eq. (H~3)) , e 2 repre- 
sents the square of the relative momentum received by 
the cc pairs per unit length of the nuclear medium. 

III. ANALYSIS OF J/i/> PRODUCTION CROSS 
SECTIONS AT FIXED TARGET ENERGIES 

In this section, we will describe the results of our 
analysis of J/ip production cross sections in hadronic 
as well as hadron-nucleus collisions. Models describing 
quarkonium production in hadronic collisions following 
color-octet or color evaporation approach generally as- 
sume that the perturbatively produced cc pairs can sub- 
sequently hadronize into observed hidden charm mesons 
(with a certain probability) provided the invariant mass 
of the pair lies in the window 2mc < Q < 2m£>. Thus 
it might be interesting to study the energy evolution of 
these sub-threshold cc pairs. In Fig. [3]we have thus plot- 
ted sub-threshold a C c[2mc, 2t71d], and total cc produc- 
tion cross section and their ratio as a function of the 
center of mass energy y/s of the colliding hadrons. At 
lower energies the two cross sections are almost indif- 
ferent indicating that all the cc pairs are predominantly 
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FIG. 3: Variation of the sub-threshold (dashed line) and total (continuous line) cc production (left panel) and their ratio (sub-threshold 
fraction) (right panel) as a function of the center of mass energy of the colliding hadrons. 



produced within the hidden charm interval. However the 
sub-threshold cross section is found to rise less steeply 
with energy compared to the total cross section. The ra- 
tio of these two cross sections, the sub-threshold fraction 
is first seen to decrease with energy then gradually sat- 
urates for large yfs. Thus higher be the energy of the 
collision, cc pairs are more likely to be produced with an 
invariant mass above the open charm production thresh- 
old. 

As pointed out earlier, at leading order, the signifi- 
cant contribution to cc production comes from the light 
quark annihilation and gluon-gluon fusions. Because of 
the larger color factors, the inclusive cross section, at a 
given energy, due to gluon fusion is always higher than 
that due to quark annihilation. However this is not the 
case if one looks at the differential cross sections of these 
two processes. To estimate their relative size, we plot 
in Fig. HI the fraction of the total differential cross sec- 
tion due to each sub-process in p-p collision for different 
values of y/s. At lower xp , gluon fusion dominates the 
cc production whereas at higher Xf, quark annihilation 
also plays a crucial role. The critical value Xp" lt , above 
which qq annihilation cross section rides over gg fusion is 
not fixed and depends on the underlying collision energy. 
Higher be the y/s, larger is the contribution due to glu- 
ons. Consequently at RHIC (y/s NN = 200 GeV) or at 
LHC (y/s NN = 5.5 TeV), J/ip production up to leading 
order can be reliably described by gg fusion alone. How- 



TABLE I: Values of fitting parameters fj/$ and qf used in 
the calculation of J/ip production cross sections 
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0.510 







1.2 GeV 


1.0 



ever, at the other end of the energy scale, in the FAIR 
energy domain, one needs to consider the light quark an- 
nihilation channel as well. This observation can be easily 
understood by looking at the behavior of parton distri- 
bution functions shown in Fig. [1] Larger xf requires a 
beam parton at large x a , which in turn means small y/s. 
In this region the gluon distribution of a proton is very 
small compared to the valence quark distribution. Con- 
sequently, the J ftp production cross section at large Xp is 
very sensitive to the sea quark distributions of the target. 

For a complete description of charmonium production, 
let us now fix the the parameters of the model, described 
in the previous section. For each of the parametric forms 
of the transition probability, combining the K-factor with 
the overall normalization constant Nj/^, we are left with 
two free parameters ap and /jm = Kj/^Nji^. The op- 
timum values of these two parameters, given in Table U 
are extracted by fitting the model results with the in- 
clusive 3/tjj production cross section data in the forward 
region (xp > 0), from pp collisions at fixed target en- 
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FIG. 4: Ratio of light quark annihilation and gluon fusion cross sections for cc production in p-p collisions to their sum as a function of 
the Feynman scaling variable xp, at different values of the center of mass energies of the collision system. MSTW 2008 LO central pdf 
set is used as an input. The solid line is for quark contribution , and the dashed line is for gluonic contribution as indicated in the figures. 



ergies. All the three parameterizations of F cE ^j/^(q 2 ) 
are independently found to fit the available data reason- 
ably well. Results are shown in Fig. [5J In addition to 
p-p reactions, data from normalized p-N collisions are 
also plotted. These nuclear target data points represent 
the forward cross section per nucleon from light targets 
(Li, Be, C). Their close agreement with the theoretical 
curves for p-p collisions indicates that for light- ions J/i/j 
production is least affected by the nuclear effects due 
to the presence of few target nucleons. A parameteriza- 
tion, called Schuler parameterization [3J for total forward 
cross section not including the branching ratio to lepton 
pairs is also available from the E672 collaboration, for 
< 31GeV, as 



a(xF > 0) = er (l — 



(14) 



where n = 12.0 ±0.9 and (Jq = 638 ±104 nb for a proton 
beam. In Fig. [5j we have also compared our results with 
this parameterization. 

Once the parameters are fixed, the model can now also 
be used to calculate the energy dependence of the cross 
section at mid-rapidity. This energy dependence has also 
been parameterized long ago by Craigie Q as: 



da -14.7m 

B— = Cexp{ j^-) 

ay V s 



(15) 



where m is the mass of the particular resonance 
state, B is the branching ratio for di-muon channel and 
C(J/V>)=40.0 nb. Fig.[B]thus shows the comparison of the 
data for the mid-rapidity J/-0 production in p-N collisions 
at fixed-target energies with the model calculations along 
with the parametric curve. 

Thus we have seen that with all three parametric forms 
of the transition probabilities, representing three possible 
physical scenarios of color neutralization, the model gives 
a reasonable description of the data for J/ip hadroproduc- 
tion at the available fixed target energies, an observation 
in line with the previous studies [l6j . Here a few words 
about the selected data sets on hadroproduction of char- 
monia may be in order. Quarkonium production cross 
sections are usually reported either as integrated over 
forward region of phase space, (xp > 0) in the center of 
mass or in the central rapidity region (y cm s = 0). The 
measurements have been performed over a time period 
spanning more than 30 years. Over such a long period, 
several experimental techniques have been used and dif- 
ferent input information were available during the time 
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FIG. 5: Total J/ip production cross section in the forward region 
(xp > 0) as a function of center of mass energy. In addition to p-p 
collisions, the plot also contains the available data points for nor- 
malized proton-nucleus collisions with light targets from different 
fixed target energies. The model parameters for each of the para- 
metric forms of the transition probability F(q 2 ) is however fixed 
from the p-p data. The cross section as obtained from Schuler 
parameterization is also included. 
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FIG. 6: Energy evolution of the mid-rapidity J/ip production 
cross section in di-muon channel for proton-nucleon collisions. The 
Craige parameterization of the differential cross section is also 
shown. Model calculations are performed following Eq. \5\ 



of measurements. Hence comparing results of different 
experiments on an equal footing requires an update of 
the published numbers on several aspects. For example, 
charmonium branching ratio has changed with time and 
the treatments of nuclear effects are not homogeneous. 
The data here are taken from [3j , where all cross section 
data are adjusted with the same di-muon branching ratio 
{Byii = 0.0597 ± 0.0025). Also nuclear target data are 
corrected to pp by assuming o p a = a pp A - 9 . 

We now examine the J/ip production for nuclear tar- 
gets. In this case the production is influenced by the 
interplay of variety of physical processes. For a thorough 
understanding of the involved mechanisms, we need to 
have accurate sets of data, spanning over large intervals 
of the incident proton beam energy and covering large 



xp and pt intervals. As far as present experimental 
situation is concerned, a wealth of high statistics J/ip 
data samples have been collected over the years in sev- 
eral fixed target experiments covering different kinematic 
range and energy domains. The HERA-B experi- 
ment at DESY has recorded J/ip production in p-C, p-W 
reactions at 920 GeV incident proton beam energy. E866 
[Io| at FNAL has studied p-Be, p-Fe and p-W collisions 
at 800 GeV, CERN-SPS experiment NA50 at 200 Q, 400 
and 450 @ GeV and NA60 QJ] at 158 and 400 GeV 
for several different nuclear targets. A unanimous fea- 
ture observed in all such measurements is that at fixed 
collision energy, charmonium production rate per target 
nucleon decreases with increasing A, the mass number 
of target nucleus. Different experimental approaches are 
in vogue to quantify the reduction. Usual way of pa- 
rameterizing the nuclear effect is to fit the A dependence 
of the production cross section with simple power law: 
c p a = Co x A a , and then to study the evolution of a 
with xf and pr- Measured values of a are found to be 
less than unity which then indicates suppression in nu- 
clear medium. Alternatively, the nuclear effects have also 
been expressed by the fitting the data in the framework of 
Glauber model having as input parameters the inelastic 
nucleon-nucleon cross section and density distributions 
for various nuclei. The model gives as an output <J^ bs , 
which represents the break-up cross section of the cc pair 
in its pre-resonance or resonance state. However, most 
often the detailed Glauber calculation, is approximated, 
in first order, by a simplified p < L > exponential pa- 
rameterization: (jpA = aoAexp(—a^ g L> p < L >). Here 
p < L > denotes the average amount of nuclear mat- 
ter crossed by the pre-formed charmonium state from 
its production point up to exiting from the nucleus. In 
NA50, at 400 GeV data set, the absorption cross sec- 
tion for J/ip has been extracted to be 4.6 ± 0.6 mb from 
Glauber fit and 4.1 ± 0.5 from p < L > fit, whereas a 
joint fit to both HI and LI data sets with the p < L > 
parameterization reproduces a value 4.3 ± 0.7 mb at 450 
GeV. By fitting the cross section ratios in the frame- 
work of Glauber model, NA60 has obtained cr'^f = 7.6 
± 0.7(stat.) ± 0.6(syst.) mb. at 158 GeV and 4.3 ± 
0.8(stat) ± 0.6(syst.) mb at 400 GeV. Recent measure- 
ments by NA60 collaboration have thus confirmed that 
the nuclear effects become more important when moving 
down the beam energy, an observation that remains valid 
when extending the comparisons to other existing sets of 
results. We should take note of the fact that both a and 
Oabs as reported by the experiments, are effective quan- 
tities, used to describe the convolution of all cold nuclear 
matter effects reducing the J/ip yield. But they do not 
allow to disentangle the different contributions (initial 
state shadowing, final state dissociation etc.) playing a 
role in this reduction. 

In the present paper, we have analyzed the J/ip pro- 
duction cross section data measured in the energy range 
from 158 to 920 GeV. The data sets from different fixed 
target experiments, that have been used in the present 
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FIG. 7: J/ip cross section ratios for p-A collisions at different energies of the incident proton beam from different fixed target experiments. 
Data are fitted with model calculations for two different forms of the transition probability representing two different physical mechanisms 
of J/ip formation from an initial cc pair. 



analysis are given in Table [TT] along with their phase space 
coverage. Among the listed experiments, NA50 presented 
their results as J/ip production cross section in di-muon 
channel for different nuclear targets as a function of the 
average nuclear path length < L >. On the other hand 
NA60 has presented their results as the ratio of J/ip pro- 
duction cross sections for a particular nucleus to that 
due to Be as a function of < L >. In case of any model 
study, one advantage to work with ratios is that we can 
get rid of the multiplicative parameters as they get can- 



celed. Hence for NA50 as well, we have fitted the model 
results with the cross section ratios rather than their ab- 
solute values. For 400 and 450 GeV we have used the 
lightest available target Be as the reference. In case of 
200 GeV, we have chosen C as our reference instead of 
H, as the nuclear suppression are found to be sensitive 
to the choice of the lightest target [2(| . Fig. [7] below il- 
lustrates our calculations fitted with the corresponding 
data sets from NA50 and NA60. On the other hand, 
at 800 (920) GeV we have fitted the ratio of 3/ij) pro- 
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TABLE II: Basic features of the experiments providing the charmonium production cross sections (or ratios) considered in the 
present study. 



Experiment 


Ijiab l^ev] 


Collision systems 


Phase space 


NA60 [12] 


158 


p-Be.ALCu, In, W, Pb, U 


0.28 < j/cms < 0.78 


NA50 [7] 


200 


p-C,Al,Cu, Ag, Pb, U 


0.0 < J/cms < 1.0 


NA50 [8] 


400 


p-Be, Al, Cu, Ag, W, Pb 


-0.425 < y cms < 0.575 


NA60 [12] 


400 


p-Be, Cu,Ag,W,Pb, U 


-0.17 < J/cms < 0.33 


NA50 [9J 


450 


p-Be,Al,Cu,Ag, W 


-0.50 < y C m S < 0.50 


E866 [10] 


800 


p-Be,W 


-0.10 < if < 0.93 


HERA-B [11] 


920 


p-C,W 


-0.34 <i F < 0.14 



TABLE III: The e 2 values obtained from fitting the data of different fixed target experiments with different energy of the 
incident proton beam. e% and e 2 P corresponds to F^ p \q 2 ) and F^°\q 2 ) respectively. The error bars correspond to fit errors. 



Experiment 


Ei ab [GeV] 


e 2 G [GeV 2 /H 


e 2 P [GcV 2 /fm] 


NA60 


158 


3.52 x 10" 1 ± 1.72 x 10" 2 


2.47 x 10" 1 ± 1.02 x 10" 2 


NA50 


200 


2.616 x 10 _1 ±8.72 x 10" 2 


1.9 -1 ± 5.67 x 10~ 2 


NA50 


400 


2.65 x lO' 1 ± 1.71 x 10~ 2 


1.86 x 10 _1 ± 1.07 x 10" 2 


NA60 


400 


2.617 x 10 _1 ± 2.22 x 10~~ 2 


1.82 x 10 _1 ± 1.37 x 10~ 2 


NA50 


450 (HI) 


1.94 x 10" 1 ± 2.46 x 10~ 2 


1.38 x 10 _1 ± 1.6 x 10" 2 


NA50 


450 (LI) 


1.79 x 10 _1 ± 3.51 x 10~ 2 


1.3 x 10 _1 ± 2.35 x 10" 2 


E866 


800 


1.83 x 10 _1 ± 3.28 x 10~ 3 


1.31 x 10 _1 ±2.18 x 10" 3 


HERA-B 


920 


1.42 x 10 _1 ± 1.19 x 10~ 2 


1.04 x 10 _1 ±8.16 x 10~ 3 



duction cross section for W to C (Be) as a function of 
xf, calculated from the parameter a and given in [26J. 
L(A), as introduced in Eq. (flU)) . is replaced by the re- 
spective value < L > for a target nucleus, as used in 
the corresponding data set at a particular fixed target 
energy. The value of e 2 (the relevant parameter in the 
model for quantifying the final state nuclear dissociation 
of J ftp) at each energy is extracted by fitting the avail- 
able data independently with the model results for both 
parameterizations of the transition probability namely 
the Gaussian distribution F^(q 2 ) (simulating the color 
singlet model) and the power law distribution F^ p \q 2 ) 
(mimicking the color octet formalism). The transition 
probability F^(q 2 ) containing the essential features of 
color evaporation model (CEM) [l8], is a special case of 
F^ p \q 2 ) for cvf = 0. It assumes that all the cc pairs 
with invariant mass less than open charm threshold has 
the same constant probability to evolve in a J/ip me- 
son. However as pointed out in the pairs just below 
the threshold should have a smaller probability to be- 
come than those far below the threshold. In other 
words, the power law parameterization F^ p \q 2 ), with 
otp > is believed represent more accurate physics for 
J/ip production. We thus refrain ourselves from further 
use of F^ c \q 2 ), while analyzing the J/tp production in 
p-A reactions. The best fit values of e 2 as obtained sep- 
arately for F^ G \q 2 ) and F^ p \q 2 ) are presented in Table 



IIIII Both £q and e 2 P show a non-negligible energy depen- 
dence with Ep/Q*, increasing with decrease in the collision 

energy. The observed energy dependence pattern of e 2 
(characterizing the final state dissociation) is thus in line 
with the earlier observations of energy dependence of the 
J/V> absorption cross section in normal nuclear matter 

It might be interesting to note that for the parameter- 
ization F^(q 2 ) in Eq. ©, a shift of q 2 to q 2 in Eq. (13]) 
for the J/ip suppression in nucleon- nucleus collisions gives 
the following relation 



(T 



;M->J/i/> 



= cxp 



2a% 



L{A) 



&NN^J/t(> 



(16) 

This relation is effectively the same as that accounted 
by the exponential p < L > parameterization. One can 
thus calculate the corresponding absorption cross section 
as <r a b 8 (m&) = (10 x e 2 )/(2 x a 2 F x p ). Using the model 
parameters at 158 GeV, and taking po ~ 0.15/ fm 3 , the 
corresponding absorption cross section comes out to be 
Cabs ~ 8.15 ± 0.4 mb. Though close, the value is still 
larger than that reported by NA60 (with a full Glauber 
model analysis). As pointed out earlier, the cr a b s as re- 
ported by different experimental groups are the effective 
quantities since they include both initial as well as final 
state effects. However in our calculation e 2 and hence 
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TABLE IV: Comparison of the extracted e 2 and the corresponding absorption cross section (a a b s ) values with free proton (p) 
and nuclear (n) parton distribution functions following EPS09 interface at different beam energies in case of Gaussian transition 
probability. The absorption cross sections extracted from the collected data from exponential p < I > fitting are also included 
in the last column. The normal nuclear matter density is taken as po = 0.15//m 3 . 



Experiment 


E lab [GeV] 


e% [GeV 2 /H 


el [GeV 2 /fm] 


°lbs [mb] 


°abs [mb] 


<bs L> [mb] 


NA60 


158 


2.96 x 10 _1 ± 1.72 x 10~ 2 


3.5 x 10 _1 ± 1.72 x 10~ 2 


6.85 ±0.4 


8.1 ±0.4 


7.36 ±0.7 


NA50 


200 


2.2 x 10" 1 ±9.4 x 10" 2 


2.62 x 10 _1 ±8.72 x 10" 2 


5.09 ± 2.08 


6.06 ±2.01 


4.68 ±3.55 


NA50 


400 


1.95 x 10 _1 ± 1.7 x 10~ 2 


2.65 x 10 _1 ± 1.7 x 10~ 2 


4.51 ±0.39 


6.13 ±0.39 


4.69 ± 0.75 


NA60 


400 


1.87 x 10 _1 ± 2.2 x 10" 2 


2.62 x 10 _1 ± 2.2 x 10" 2 


4.3 ±0.51 


6.06 ±0.51 


4.44 ± 1.02 


NA50 


450 (HI) 


1.5 x 10 _1 ±2.6 x 10~ 2 


1.94 x 10 _1 ± 2.46 x 10~ 2 


3.47 ±0.6 


4.49 ±0.57 


3.8 ± 1.1 


NA50 


450 (LI) 


1.34 x 10 _1 ±3.8 x 10" 2 


1.79 x 10 _1 ±3.51 x 10" 2 


3.1 ±0.88 


4.14 ±0.81 


3.71 ± 1.51 
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FIG. 8: Energy evolution of the shadowing factor S^ b at mid- 
rapidity, following EPS09 parameterization. The individual contri- 
bution due to light quark annihilation (open triangles) and gluon 
fusion (open crosses) are shown separately along with their joint 
contribution (open circles). 



the corresponding cr a bs characterizes only the size of the 
final state dissociation. We can calculate the effective e 2 
(forf( G >(g 2 )) at different energies, following our adopted 
model, if we neglect the initial state effect in our cal- 
culation and use free proton pdfs without any nuclear 
modification for calculation of cc pair production cross 
section. Table IIVI shows a comparison between the best 
fitted e 2 values at different SPS energies with and with- 
out the nuclear modification of the parton distribution 
functions and along with the corresponding <7 a b s - In ad- 
dition, the effective absorption cross sections extracted 
from the exponential fitting of the experimentally mea- 
sured data points are also shown. 

It appears from the Table IIVI that if we incorporate, 
as initial state effect, the nuclear modification of the par- 
ton densities following EPS09 interface, we need a higher 
value of e 2 than that of a free proton pdf to reproduce the 
corresponding data, in the SPS energy regime. This can 
be understood by studying the influence of the nuclear 
modification of the parton distribution functions to the 
J/ip production. For this purpose let us define the shad- 
owing factor {Sji^) for a particular nucleus as the ratio 
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FIG. 9: Variation of the shadowing factor as function of < 

L > at NA60 for two different energies Ei, = 158 GeV (upper panel) 
and Ei, = 400 GeV (lower panel) of the incident proton beam. The 
contribution of the light quark annihilation and gluon fusion are 
separately shown along their sum. The inclusive production cross 
sections are obtained by integrating the differential cross sections 
over the respective rapidity windows corresponding to the two beam 
energies. The Gaussian form of transition probability is used for 
evolution of cc pairs into J /if) mesons. 



between the i/ip production cross section (per nucleon- 
nucleon) collision in p-A and p-p: 
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Thus SKf, less than unity indicates shadowing whereas 



S P u^ larger than unity signifies anti-shadowing. Since 
the pdfs enter as an input to the calculation of pertur- 
bative cc pair production cross section, hence in FigJSl 
we plot S^? b at mid-rapidity, as a function of center of 
mass energy, in the domain of different fixed target ex- 
periments. In the energy regime probed by the SPS ex- 
periments, charm production at mid-rapidity explores x 
(fraction of nucleon momentum carried by partons) val- 
ues corresponding to the anti-shadowing region, where 
parton densities in the nuclei are enhanced with respect 
to those of free protons. In the absence of any final state 
interaction this would certainly result in enhancement of 
3 ftp production cross section per nucleon in p-Pb colli- 
sions compared to p-p collisions. However in the lower 
energy regime, relevant for FAIR experiments, EPS09 
parameterization shows a shadowing effect resulting in 
the reduction of the produced cc pairs in p-Pb collisions 
compared to the p-p collisions at a given energy. 

To make our study more quantitative let us now cal- 
culate for different nuclei as a function of < L > 
at NA60 beam energies, corresponding to their respec- 
tive kinematic interval. The results are depicted in Fig. 
(J9j> . As it appears from the figure, at both the energies 
we obtain an over all anti-shadowing for all the nuclei, 
resulting in an increment in the 3 ftp production cross 
section compared to the p-p reactions in absence of any 
final state dissociation. At 158 GeV, the qq fraction of 
the cross section shows a very weak shadowing for heav- 
ier nuclei, which is more than counterbalanced by the 
stronger contribution from gluon fusion. Hence if nuclear 
modification of parton densities is taken into account, 
one needs a more severe amount of final state multiple 



scattering, than with a free proton pdf, to reproduce the 
measured data collected at SPS, which actually show a 
suppression in yield as a function of < L >. Naturally 
the absorption cross section extracted from the corre- 
sponding e 2 will also be larger if one takes into account 
the nuclear pdfs following EPS09 interface compared to 
the free proton (as well as the experimentally measured) 
case. The observation is in agreement with the earlier 
studies [H, • The results are sensitive to the adopted 
parameterization of the nuclear parton densities. For ex- 
ample, the NA60 collaboration has observed that if the 
initial state effects are evaluated with EKS98 [3(| pa- 
rameterization, a larger absorption cross section is re- 
quired to fit the measured data (c^f (158 GeV) = 9.3 
FIG. 10: Variation of the shadowing factor S*^ as function of < ±07±0.7 mb and a J J^ (400 GeV) = 6.0 ± 0.9 ± 0.7 mb). 

L > for a 30 GeV incident proton beam. Such measurements will be Slightly higher a'tf has been obtained Using EPS08 [3l| 
performed in the SIS-100 phase of the upcoming CBM experiment , . ,. TT , n . i, i j. £ 

at FAIR parameterization. However the situation mignt be uit- 

ferent at FAIR as seen from Fig. [TO] The inclusive pro- 
duction cross sections are calculated over a rapidity in- 
terval —0.5 < ycms < 0.5. Incorporation of nuclear pdfs 
following EPS09 parameterization, results in reduction 
of the cc and hence 3 ftp production in nuclear targets 
compared to the protonic collisions. These observations 
can be understood by looking at the nature of EPS09 
parameterization in Fig. [5] The rapidity coverage of 
NA60 at 158 GeV is 0.28 < y cms < 0.78. This roughly 
corresponds to a Xb{= {mj/^/y/s)exp(—y cms )) regime, 
~ 0.08 < x\j < 0.13, where we observe anti-shadowing. 
However the rapidity interval at 30 GeV beam energy 
considered in the present calculation explores a Xb region, 
~ 0.24 < Xb < 0.66. In this region parton population in 
a nucleon bound inside a nucleus is depleted compared 
to a free nucleon. Thus, within the present theoretical 
scenario, a smaller dissociation cross section, compared 
to the free proton case might be required to describe the 
3/ip suppression expected to occur at FAIR. However a 
conclusive picture can be drawn only after the measure- 
ments are performed. 

Apart from the inclusive production cross sections, 
NA60 has also measured the 3/tp rapidity distributions 
[lH for several target nuclei at 158 GeV as well as 400 
GeV. Within the experimental coverage of NA60 di- 
muon spectrometer, the rapidity distributions are found 
to be well reproduced by Gaussian functions. Prelimi- 
nary analysis reports that at 158 GeV, the distributions 
corresponding to the different p-A collisions can simul- 
taneously be described, by a single Gaussian function, 
having mean (fi y ) = 0.05 ± 0.05 and sigma (a) = 0.51 ± 
0.02. At 400 GeV the corresponding fit parameters are 
mean {n y ) ~ -0.2 (guided by the previous NA50 mea- 
surements at the same energy) and (a) = 0.81 ± 0.03. 
However the data for rapidity distribution are presented 
in arbitrary units. Cross section values have not yet been 
published. We thus compare the shape of rapidity dis- 
tributions at both energies, as obtained from our cal- 
culations with the experimentally measured ones. The 
rapidity distributions for different collision systems are 



12 




0.5 
y 

' cms 




-0.5 0.5 



' cms 

FIG. 11: 3/ip rapidity distributions in p-A collisions for Be and 
Pb targets at 158 GeV (top) and 400 GeV (bottom) as obtained in 
our calculations. Two different parametric forms of the transition 
probabilities produce almost similar results. 

simultaneously fitted by a single three parameter Gaus- 
sian function. For an unbiased comparison we have also 
performed the Gaussian fitting of the experimental data 
points. Results are given in Table [V] At both the ener- 
gies, the widths of the rapidity distribution are found to 
be higher in our model calculations in comparison with 
the measured values. Moreover our calculations can not 
reproduce the negative mean value of the rapidity distri- 
bution at 400 GeV. For illustration, we have shown the 
rapidity distributions for Be and Pb targets, following 
our model calculations in Fig. QTJ 



IV. 3/i> PRODUCTION AT FAIR 

Our ultimate goal is to estimate the 3/ip yield in p-p 
and p-A collisions at energy values planned to be ex- 
plored at FAIR. In the start version of FAIR (SIS-100), 
proton beams up to 30 GeV will be accelerated. In the 
final version high intensity proton beams with energy up 
to 45 GeV will be available. In Fig. [12l we present the 
J/V> rapidity distribution in protonic collisions with 15, 
30 and 45 GeV proton beams, for all three different func- 
tional forms of the transition probability. As expected, 
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FIG. 12: Rapidity distribution of ]/ip mesons produced in pp 
collisions in di-muon channel, for 15 GeV (upper panel), 30 GeV 
(middle panel) and 45 GeV (bottom panel) incident proton beams. 
Such measurements will be performed in the SIS-100 and SIS-300 
rings of the upcoming CBM experiment at FAIR. 

the distributions are symmetric about the center of mass 
rapidity y cms = 0, and the width of the distribution de- 
creases with the decrease in the beam energy. The data 
for inclusive J/?/> production cross section, in this energy 
regime, are well accounted by all three different cases of 
formation mechanisms as seen in the previous section. 
However the differential cross sections take different val- 
ues for different form of F cS ^j/ 1 p(q 2 ). The yield is found 
to be minimum for F( c )(q 2 ) and maximum for F^ p \q 2 ) 
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TABLE V: Parameters of the Gaussian functions used to fit the J/ip rapidity distributions in p-A collisions at 158 GeV and 
400 GeV. The superscripts (P) and (G) respectively denote the power law and Gaussian parameterizations. 



E Lab [GeV] 


(Expt) 




(G) 
My 


(Expt) 
°V 


4 P) 


4 G) 


158 


0.044 ±0.1 


0.02 ± 0.087 


0.02 ±0.086 


0.55 ±0.06 


0.66 ±0.09 


0.065 ± 0.096 


400 


-0.27 ±0.16 


9.17 x 10^ 4 ±0.037 


-1.93 x 10" 4 ±0.034 


0.86 ±0.2 


0.97 ±0.07 


0.92 ±0.059 
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FIG. 13: Variation of e 2 with energy of the incident proton beam 
for a) F< p )(q 2 ) (top panel) and b) F^ G \q 2 ) (bottom panel). The 
points on the plots are taken from IIIII In both the cases the 
energy dependence has been parameterized following exponential 
functions. 



for all three beam energies, with the relative difference 
decreasing with increasing beam energy. 

Next we present results for charmonium production in 
p-A collisions at FAIR. For this purpose, we have param- 
eterized the observed dependences of e 2 on the energy of 
the incident proton beam for both F^(q 2 ) and F<- P \q 2 ) 
using exponential functions. The parametric equations 
£ 2 (ELab) are shown in Fig. [13l These parametric rela- 
tions can now be extrapolated to derive the level of sup- 
pression at energies relevant to the CBM experiment at 
FAIR. In the Fig. [14] we have predicted the normal nu- 
clear suppression in p-A collisions, for E^ab =15, 30 and 
45 GeV. Following the same approach as NA60, we have 
calculated J/ip cross section ratios (with Be as the refer- 



ence) as a function of < L > for different target nuclei. 
In case of F^- GS> (q 2 ), we can also find the corresponding 
absorption cross section as might be extracted from the 
data that will be collected at FAIR. Results are given 
in Table IVII However the experimentally measured val- 
ues can be larger than those reported here due to parton 
shadowing as explained in the last section. In addition 
to the inclusive production cross sections, we also cal- 
culate the J/ip rapidity distributions for different target 
nuclei. Fig. [T5] illustrates the J/ip differential production 
cross section at these three energies for three different 
p-A systems. The obtained distributions can be fitted 
with the Gaussian functions. The corresponding best fit 
parameters are given in the Table [VlTJ 



TABLE VI: J/ip absorption cross section (a a bs) in cold nuclear 
matter at FAIR energies. EPS09 parameterization has been 
used as an input for calculation of J/ip production and thus 
the dabs characterizes only final state dissociation. 



Ei ab [GeV] 


a abs [mb] 


15 


10.4 ± 1.79 


30 


10.1 ± 1.77 


45 


9.87 ± 1.76 



At SPS and higher energies, the calculated inclusive 
J/ip yields employing both the forms of transition prob- 
abilities, are found to match the measured data for dif- 
ferent p-A systems, within error bars. Results for dif- 
ferential distributions are also in close proximity for the 
two cases. In other none of the inclusive or differen- 
tial measurements help us to select either of the exist- 
ing schemes of color neutralization. However the pic- 
ture appears to be different at lower energies. The two 
possible hadronization schemes are found to generate 
different amount of suppressions. The Gaussian distri- 
bution gives higher suppression compared to the power 
law distribution. As the energy of the incident proton 
beam decreases, the difference between the J/ip yields 
as predicted by two models increases. This might be 
attributed to the difference between the physical mecha- 
nisms incorporated in F^ G \q 2 ) and F^(q 2 ). The phase 
space between the thresholds of producing a cc pair and 
open charm meson is fairly large (~ 5 GeV 2 ). Conse- 
quently semi-hard gluons can be radiated during the for- 
mation of cc pairs with large invariant mass. The ra- 
diation reduces the invariant mass of the pair and thus 
strongly enhances the probability for the pair to form 
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FIG. 14: J/ip production cross section ratios as a function of < L > for different target nucleus in the FAIR energy domain, as predicted 
by our calculations. We have used both the Gaussian F^ G \q 2 ) as well as power law F^ p '(q 2 ) forms of the transition probabilities for 
calculation of J/tp yield. Lower be the energy of the incident proton beam the two different forms representing two different physical 
mechanisms of J/ip formation is seen to generate different amount of suppression. 



a J/ip in case of F^(q 2 ). However, for the Gaus- 
sian probability F^^q 2 ), the related physical process 
is not assisted by any gluon radiation during the evolu- 
tion. As a consequence, the J/ip formation probability 
is highly suppressed for the cc pairs produced with large 
invariant mass. This results in a difference between the 
J/ip production cross sections following two hadroniza- 
tion schemes, which becomes appreciable at lower beam 
energies. This can be understood by looking at the in- 
variant mass distribution of the produced cc pairs at dif- 
ferent beam energies. Fig. [T5] shows the distribution of 
the cc pairs produced at the initial stage (from gluon 
fusion and light quark annihilation) at the two energies 
E Lab = 30 GeV and 400 GeV. The shape of the distri- 
butions are quite different in the two cases. At 400 GeV 
the produced pairs are almost uniformly distributed over 
the allowed phase space interval. However at 30 GeV, 
the production of large mass cc pairs is highly depleted 
due to the lack of center of mass energy. Each of these 
input distributions are convoluted with the two differ- 
ent functional forms, F^ p \q 2 ) and F^(q 2 ), to obtain 
the final state distributions, as presented in Fig. [T7] At 
lower energy, the separation between the two distribu- 
tions, corresponding to two different physical processes 



of hadronization, becomes higher. This ultimately re- 
sults in the different production cross sections for two 
different processes. At FAIR very high statistics J/ip 
data will be collected over a fairly large period of time. 
Hence the measurements are expected to be least suf- 
fered by the statistical uncertainties. Thus in addition to 
the accurate normalization of the heavy-ion data, precise 
measurements of inclusive as well as differential J/ip pro- 
duction cross sections in p-p and p-A collisions at FAIR 
might also help to distinguish the predictions for different 
scenarios and thereby pinning down the correct mecha- 
nism for color-neutralization. 



V. SUMMARY 

The study of charmonium production in high energy 
collisions is highly advantageous. Because of the large 
charm quark mass (mc » ^qcd), and the limitation 
to be produced at pair (to conserve flavor in strong in- 
teraction) , charmed hadrons provide us with a new scale 
to probe QCD. In the present article, we have made a 
systematic study of J/ip production in proton induced 
reactions, in the beam energy range from 158 GeV to 



15 




FIG. 15: The differential J/-0 production cross section per nucleon as a function of the center of mass rapidity y cm s for different target 
nucleus (Be, In and Pb) in the FAIR energy domain, at 'EiLab = 15 GeV (top), 30 GeV (middle), 45 GeV (bottom), as predicted by our 
calculations. Both the Gaussian F^ G \q 2 ) (markers) as well as power law F^ p ~>(q 2 ) (lines) forms of the transition probabilities are used 
for calculation of J/tp yield. 



TABLE VII: Parameters of the Gaussian functions used to fit the J/ip rapidity distributions in p-A collisions for incident proton 
beam energies 15 GeV, 30 GeV and 45 GeV. The superscripts (P) and (G) respectively denote the power law and Gaussian 
parameterizations. 



E Lab [GeV] 


(p) 
fJ-y 


(G) 
M 


Uy 


Uy 


15 


1.34xl0~ 2 


1.28xl0~ 2 


2.11xl0 _1 


2.04xl0 _1 


30 


1.81xl0~ 2 


1.82xl0~ 2 


3.57xl0 _1 


3.51xl0 _1 


45 


1.78xl0~ 2 


1.89xl0~ 2 


4.52xl0 _1 


4.41xl0~ 1 



920 GeV, available from different fixed target experi- 
ments. The model employed for our study is based on the 
QCD factorization of inclusive processes. It is general 
enough to incorporate various existing physical schemes 
of color neutralization through parametric functions. In 
case of p-A collisions, two major cold nuclear matter 
effects namely the modification of parton densities in- 
side the target nucleons at the initial stage and the fi- 
nal state multiple scattering of the nascent cc pairs in- 
side the target are taken into consideration. Interplay of 
these two effects on the measured J/ip production cross 
section has been investigated. However a more rigorous 
theoretical framework should also include other possible 



nuclear matter effects like parton energy loss at the ini- 
tial and final stage, that may influence the production. 
The model gives a satisfactory description of the avail- 
able data on J/ip production cross sections in nucleon- 
nucleon and nucleon-nucleus collisions. However such 
measurements are not found suitable to discard any of 
the existing physical mechanisms of color neutralization, 
when compared to the data. The final state interaction 
shows a non-trivial dependence on the energy of the in- 
cident proton beam. Lower be the energy of the pro- 
jectile proton, more prominent is the effect of the nu- 
clear medium on J/ip production. This observation is in 
agreement with the existing theoretical and experimen- 
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FIG. 16: The invariant mass distribution of the initially produced cc pairs at E^ a (, = 30 GeV and 'EiLab = 400 GeV. 




FIG. 17: The invariant mass distribution of the cc pairs at ^Lab = 30 GeV and E^ a f, = 400 GeV, at the final stage after they are 
convoluted with the transition probabilities F( p )(g 2 ) and F^ G \q 2 ) for J/ip formation. 



tal observations and has a direct impact to the heavy-ion 
collisions. At lower beam energies, the colliding nuclei 
are less Lorentz contracted and thus results in larger col- 
lision time (t co u = 2Ra/j) as well as formation time 
of the partonic medium. Thus the J/"0 mesons, during 
its evolution, are more likely to encounter the (primary) 
nuclear medium rather than the (secondary) de-confined 
phase. Results of our present investigations are utilized 
to make predictions about the level of J ftp production 
in p-p and p-A collisions at energies, close to the kine- 
matic threshold. The future FAIR experiments will take 
data in this energy domain. An accurate measurement 
of charmonium production production in p-p and p-A 



collisions performed in the same kinematic conditions 
of A-A collisions, appears to be an integral part of the 
charmonium program at FAIR. A proper calibration of 
the conventional cold nuclear matter effects will help us 
to identify the unconventional anomalous suppression in 
nucleus-nucleus collisions that might result from the on- 
set of color de-confinemcnt in the produced system at 
high baryon densities. A precision measurement of the 
charmonium production cross sections for different nu- 
clear targets, may also help us to shed light on the much 
debated issue of color neutralization of a cc pair leading 
to the formation of hidden charm mesons. 
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